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o 

' Spin and chirality orderings of the three-dimensional Heisenberg spin glass under magnetic fields 

' are studied by large-scale equilibrium Monte Carlo simulations. It is found that the chiral-glass 

transition and the chiral-glass ordered state, which are essentially of the same character as their 
O ' zero-field counterparts, occur under magnetic fields. The chiral-glass ordered state exhibits a one- 

I step-like peculiar replica-symmetry breaking in the chiral sector, while it does not accompany the 

spin-glass order perpendicular to the applied field. Critical perperties of the chiral-glass transition 
are different from those of the standard Ising spin glass. Magnetic phase diagram of the model is 
constructed, which reveals that the chiral-glass state is quite robust against magnetic fields. The 
chiral-glass transition line has a character of the Gabay- Toulouse line of the mean-field model, yet its 
' physical origin being entirely different. These numerical results are discussed in light of the recently 

C , developed spin-chirality decoupling-recoupling scenario. Implication to experimental phase diagram 

' is also discussed. 
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I. INTRODUCTION 



In the studies of spin glasses, much effort has been devoted either exprimentally or theoretically to the properties 
' O under magnetic fields. Unfortunately, our understanding of them still has remained unsatisfactory On theoretical 
C ' side, most of the numerical studies have focused on the properties of the simple Ising model, especially the three- 
dimensional (3D) Edwards-Anderson (EA) model. While the existence of a true thermodynamic spin-glass (SG) 
transition has been established for this model in zero field, the question of its existence or nonexistence in magnetic 
fields has remained unsettled. This question is closely related to the hotly debated issue of whether the ordered state 
of the 3D Ising SG in zero field exhibits a replica-symmetry breaking (RSB) or not. 

If one tries to understand real experimental SG ordering, one has to remember that many of real SG materials 
are more or less Heisenberg-like rather than Ising, in the sense that the random magnetic anisotropy is considerably 

■ weaker than the isotropic exchange interaction For example, in widely studied canonical spin glasses, i.e., dilute 
metallic alloys such as AuFe, AgMn and CuMn, random magnetic anisotropy originated from the Dzyaloshinski-Moriya 

T-H interaction or the dipolar interaction is often one or two magnitudes weaker than the isotropic RKKY interaction. 
1 Numerical simulations have indicated that the isotropic 3D Heisenberg SG with finite-range interaction does not 

■ exhibit the conventional SG order at finite temperature in zero field (However, see also Ref. |^.) Since applied 
' fields generally tend to suppress the SG ordering, a true thermodynamic SG transition is even more unlikely under 

^ magnetic fields in case of the 3D Heisenberg SG. 

Ch ' Experimentally, however, a rather sharp transition-like behavior has been observed under magnetic fields in typical 
Heisenberg-like SG magnets, although it is not completely clear whether the observed anomaly corresponds to a true 
(-H thermodynamic transition [l| JlO| , pd[ | . The situation is in contrast to the zero-field case where the existence of a true 
Q thermodynamic SG transition has been established experimentally [0. Set aside the question of the strict nature of 
O ' the SG "transition", it is experimentally observed that a weak applied field suppresses the zero-field SG transition 
^ ' temperature rather quickly. For higher fields, the SG "transition" becomes much more robust to fields, where the 
"transition temperature" shows much less field dependence p|, p^ , pl| . Such behaviors of the SG transition temperature 
under magnetic fields Tg(H) were often interpreted in terms of the mean-field model Indeed, the mean-field 

Sherrington-Kirkpatrick (SK) model with an infinite-range Heisenberg exchange interaction with weak random 
magnetic anisotropy exhibits a transition line similar to the experimental one JTst , i.e., the so-called de Almeida- 
Thouless (AT) line oc {Tg{Q) - Tg{H)f/^ in weak-field regime where the anisotropy is important, and the 
Gabay-Toulouse (GT) line |l| i? oc (Tg(0) - Tg{H))^/^ in strong-field regime where the anisotropy is unimportant. 
Nevertheless, if one notes that the true finite-temperature transition under magnetic fields, though possible in the 
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infinite-range SK model, is unlikely to occur in a more realistic finite-range Heisenberg model, an apparent success of 
the mean-field model in explaining the experimental phase diagram should be taken with strong reservation. Thus, the 
question of the true nature of the experimentally observed SG "transition" under magnetic fields remains unsolved. 

Recently, one of the present authors has proposed a scenario, the spin-chirality decoupling-recoupling scenario, 
aimed at explaining some of the puzzles concerning the experimentally observed SG transition In this scenario, 
chirality, which is a multispin variable representing the sense or the handedness of local noncoplanar spin structures 
induced by spin frustration, plays an essential role. As illustrated in Fig.|l|, locally noncoplanar spin structures 
inherent to the SG ordered state sustain two energetically degenerate "chiral" states, "right-handed" and "left- 
handed" states, characterized by mutually opposite signs of the "chiralities" . Here, one may define the local chirality 
by three neighboring Heisenberg spins Si, S2 and ^3 by, 

X-5i- (52X^3) . (1) 

This type of chirality is called "scalar chirality", in distinction with "vector chirality" defined as a vector product of 
two neighboring Heisenberg spins, S^i x 5^2 Note that the chirality defined by Eq.(|l|) is a pseudoscalar in the 

sense that it is invariant under global SO{3) spin rotations but changes its sign under Z2 spin reflections (or inversions 
which can be viewed as a combination of reflections and rotations). 




For a fully isotropic Heisenberg SG, in particular, the chirality scenario of Ref. y| claims the occurrence of a 
novel chiral- glass ordered state in which only the chirality exhibits a glassy long-range order (LRO) while the spin 
remains paramagnetic. At the chiral-glass transition, among the global symmetries of the Hamiltonian, 0(3) = 
Z2 X 50(3), only the Z2 spin reflection (inversion) symmetry is broken spontaneously with keeping the SO{3) spin 
rotation symmetry preserved. Note that this picture entails the spin-chirality (or 5*0(3) — Z2) decoupling on long 
length and time scales: Namely, although the chirality is not independent of the spin on microscopic length scale, it 
eventually exhibits a long-distance behavior entirely different from the spin. Such a chiral-glass transition without 
the conventional spin-glass order was indeed observed in recent equilibrium and off-equilibrium Monte Carlo (MC) 
simulations in zero field performed by Hukushima and one of the authors (H.K.) 0,^. It was also found there that 
the critical properties associated with the chiral-glass transition were different from those of the Ising SG, and that 
the chiral-glass ordered state exhibited a one-step- like peculiar RSB. 

In the chirality scenario of Ref. experimental SG transition in real Heisneberg-like SG magnets is regarded 
essentially as a chiral-glass transition "revealed" via the random magnetic anisotropy. Weak but finite random 
magnetic anisotropy inherent to real magnets "recouples" the spin to the chirality, and the chiral-glass transition 
shows up as an experimentally observable spm-glass transition in real Heisneberg-like SG magnets. An interesting 
outcome of this picture is that the experimental SG transition is dictated by the chiral-glass transition of the fully 
isotropic system, not by the spin-glass transition of the fully isotropic system, which has been separated from the 
chiral one. 

Very recently, the present authors discussed some of the possible consequences of the chirality scenario of Ref. Q| 
on the finite-field properties of the fully isotropic 3D Heisenberg SG [0. It was argued there that the chiral-glass 
transition, essentially of the same character as the zero-field one, occurred also in finite fields. In the weak field regime, 
the transition line was predicted to behave as 
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Tcg(O) - Tcg{H) =cH^ + dH^ + --- 



(2) 



where c and d are constants. Generally, the coefRcient c could be either positive or negative. An interesting observation 
here is that the chiral-glass transition line (Q) aparrently has a form similar to the GT line of the mean-field model. 
We emphasize, however, that their physical origin is entirely different. The quadratic dependence of the chiral-glass 
transition line is simply of regular origin, whereas that of the GT-line in the SK model cannot be regarded so. 

In the present paper, we report on our results of large-scale Monte Carlo simulations on the 3D isotropic Heisenberg 
SG, performed with the aim to reexamine the SG ordering in magnetic fields in light of the chirality scenario. In 
particular, by means of extensive numerical simulations, we wish to clarify in detail how the spin and the chirality 
order in applied fields. Part of the MC results have been reported in Ref. p^ . 

The present paper is organized as follows. In §||, we introduce our model and explain some of the details of our 
numerical method. Various physical quantities calculated in our MC simulations are defined in §111. The resu lts of 
MC s imula tions are presented in §|^. The results for the chirality- and spin-related quantities are presented in § [V A 



and §[VB, respectively. It is found that the chiral-glass transition, essentially of the same character as the zero-field 
one, occurs under magnetic fields. The chiral-glass ordered state exhibits a one-step-like peculiar replica-symmetry 
breaking in the chiral sector, while it does not accompany the s pin-g lass order perpendicular to the applied field. 
Critical properties of the chiral-glass transition are analyzed in § |IV C| . The analysis suggests that the universality 
class of both the zero-field and finite- field chiral-glass transitions might be common, which, however, differs from 
that of the standard 3D Ising SG. In §[VD, we construct a magnetic phase diagram of the model. The chiral-glass 
ordered state remains quite robust against magnetic fields, while the chiral-glass transition line in applied fields has 
a character of the GT line of the mean-field model. Section ^ is devoted to summary and discussion. Our numerical 
results are discussed in terms of the recent experimental result on canonical SG. 
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II. THE MODEL AND THE METHOD 



In this section, we introduce our model and explain some of the details of our numerical method. The model we 
consider is the isotropic classical Heisenberg model on a 3D simple cubic lattice defined by the Hamiltonian, 



where Si — {Sf,Sf, Sf) is a three-component unit vector, and H is the intensity of magnetic field applied along the z 
direction. The nearest-neighbor coupling Jy is assumed to take either the value J or — J with equal probability (± J 
distribution). 

We perform equilibrium MC simulations on this model. Simulations are performed for a variety of fields H/ J — 0.05, 
0.1, 0.5, 2.0, 3.0, and 5.0, while most extensive calculations are performed for H/ J = 0.1 and 0.5. The lattices studied 
are simple-cubic lattices with N = sites with L = 6, 8, 10, 12 and 16 with periodic boundary conditions. Sample 
average is taken over 128-1400 independent bond realizations, depending on the system size L and the field intensity 
H. Limited amount of data is also taken for L = 20 in some cases (30 samples only) to check the size dependence of 
some physical quantities. To facilitate efficient thermalization, we combine the standard heat-bath method with the 
temperature-exchange technique |l^ . Care is taken to be sure that the system is fully equilibrated. Equilibration is 
checked by the following procedures: First, we monitor the system to travel back and forth many times during the 
the temperature-exchange process (typically more than 10 times) between the maximum and minimum temperature 
points, and at the same time check that the relaxation due to the standard heat-bath updating is reasonably fast at 
the highest temperature, whose relaxation time is of order 10^ Monte Carlo steps per spin (MCS). This guarantees 
that different parts of the phase space are sampled in each "cycle" of the temperature-exchange run. Second, we 
check the stability of the results against at least three times longer runs for a subset of samples. Error bars of physical 
quantities are estimated by the sample-to-sample statistical fluctuation over bond realizations. Further details of our 
Monte Carlo simulations are given in Table |. 




(3) 
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TABLE I. Details of our MC simulations. H/ J represents the magnetic-field intensity, L the lattice size, A's the total 
number of samples, iVr the total number of temperature points used in the temperature-exchange run, Tmax/ J and Tmin/ J the 
maximum and minimum temperatures in the temperature-exchange run. 



H/J 


L 




iVT 








6 


800 


26 


0.40 


0.085 


0.05 


8 


600 


20 


0.40 


0.15 




10 


184 


34 


0.40 


0.155 




12 


56 


34 


0.40 


0.155 




6 


1400 


26 


0.40 


0.085 




8 


1300 


30 


0.40 


0.095 


0.1 


10 


1066 


46 


0.40 


0.0975 




12 


680 


56 


0.40 


0.0975 




16 


180 


50 


0.32 


0.12 




(20) 


(30) 


(46) 


(0.30) 


(0.15) 




6 


800 


26 


0.40 


0.085 




8 


800 


20 


0.40 


0.15 


0.5 


10 


552 


34 


0.40 


0.155 




12 


456 


34 


0.40 


0.155 




16 


128 


36 


0.32 


0.153 




(20) 


(26) 


(46) 


(0.30) 


(0.15) 




6 


800 


26 


0.40 


0.085 


2.0 


8 


800 


20 


0.40 


0.15 




10 


92 


34 


0.40 


0.155 




6 


800 


26 


0.40 


0.085 


3.0 


8 


800 


20 


0.40 


0.15 




10 


64 


34 


0.40 


0.155 




6 


800 


26 


0.40 


0.085 


5.0 


8 


800 


20 


0.40 


0.15 




10 


92 


34 


0.40 


0.155 
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III. PHYSICAL QUANTITIES 



In this section, we define various physical quantities calculated in our simulations below. 



A. Chirality-related quantities 

Let us begin with the definition of the chirality. We define the local chirality at the i-th site and in the /i-th 
direction, Xifn for three neighboring Heisenberg spins by the scalar 

X^, = S,^e,■iS,xS^^^^ (4) 

where {fi = x, y, z) denotes a unit vector along the /x-th axis. By this definition, there are in total 3A^ local chiral 
variables in the system. 

The mean local amplitude of the chirality, may be defined by 

1 ^ 

-2 



E Kxi>] , (5) 



where (• • •) represents the thermal average and [• • •] the average over the bond disorder. This quantity vanishes for 
coplanar spin structures, and its magnitude tells us the extent of the noncoplanarity of the local spin structures. 

By considering two independent systems ("replicas") described by the same Hamiltonian (^, one can define an 
overlap of the chiral variable via the relation, 

'^x=^E E , (6) 

1=1 fj,=x,y,z 

where Xil^ xfj represent the chiral variables of the replicas 1 and 2, respectively. In our simulations, we prepare 
the two replicas 1 and 2 by running two independent sequences of systems in parallel with diff'erent spin initial 
conditions and different sequences of random numbers. In terms of this chiral overlap g^, the chiral-glass order 
parameter may be defined by 

while the associated chiral-glass susceptibility may be defined by 

X^ = 3N[{ql)] . (8) 

Unlike the spin variable, the local magnitude of the chirality is temperature dependent somewhat. In order to take 
account of this short-range order effect, we also consider the reduced chiral-glass order parameter and the reduced 

(2) _ 

chiral-glass susceptibility Xx by dividing and Xx appropriate powers of x-i 

^(2) ^ 

X X 

The Binder ratio of the chirality is defined by 

5x = :^(3-j^) . (10) 
One may also define the distribution function of the chiral overlap by 

Px{i'x)-mq'x-<ix))] ■ (11) 

In order to study the equilibrium dynamics of the model, we also compute the autocorrelation function of the 
chirality defined by 
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1 ^ 

^xW = ^E E [(x.p(io)x^p(i + M>] ■ (12) 

i—l ^—x.y,z 

where the "time" t is measured in units of MCS. In computing simulation is performed according to the standard 
heat-bath updating without the temperature-exchange procedure, while the starting spin configuration at t = to is 
taken from the equilibrium spin configurations generated in our temperature-exchange MC runs. 

We also calculate the so-called G and A parameters for the chirality, recently discussed in the literature [p^[-p5[, 
defined by, 



- r/„2\i2 FTlfn' (13) 





- 




- U)] 




[{4)] 




J2 



Ax = " r/l2M2 " • (14) 

These and parameters are closely related to the sample-to-sample fluctuation of the chiral order parameter. 
The A parameter is known to be an indicator of the non-self-averagingness of the order parameter, i.e., it vanishes 
in the state where the order parameter is self-averaging and takes a nonzero value otherwise By contrast, the 
G parameter could take a nonzero value even in a self-averaging ordered state, and hence, cannot be used as an 
unambiguous indicator of the non-self-averagingness ||20[| . However, since in the thermodynamic limit it vanishes in 
the high-temperature phase and takes a nonzero value in the ordered state, it can still be used as an indicator of a 
phase transition. 

B. Spin- related quantities 

As in the case of the chirality, it is convenient to define an overlap variable for the Heisenberg spin. In this case, 
the overlap might naturally be defined as a tensor variable q^^ between the /x and v components (/it, v=x, y, z) of the 
Heisenberg spin. 



1 ^ 

'ii^^ = J^^^'w'^il^ ' (A' = 2;,y,z), (15) 

i=l 

where s'f'^ and sf''' are the i-th Heisneberg spins of the replicas 1 and 2, respectively. In terms of these tensor 
overlaps, the "longitudinal" (parallel to the applied field) and "transverse" (perpendicular to the applied field) SG 
order parameters may be defined by 

'Zi'^ - ll = <ll. (16) 

= [(4)], 4= E ll'^ ■ (17) 
The associated longitudinal and transverse Binder ratios are defined by 



2 V I(9£>1' 



2 V Mil)? 

Here, gi^ and are normalized so that, in the thermodynamic limit, they vanish in the high-temperature phase and 
gives unity in the nondegenrate ordered state. 

Since an odd quantity like ((/l) does not vanish in applied fields, one can also define the "connected Binder ratio" 
for the longitudinal component |24[ , 
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,i 'A-j-^Ma,, ,20, 



2 V l«9I, - M)')? 

In applied fields, might well behave differently from gi^. 

The spin-overlap distribution function is generally defined in the tensor space. In the following, we pay particular 
attention to its transverse (XY) part. The relevant transverse overlap originally has 2^ — 4 independent components. 
For the convenience of illustration, we follow Ref. [p6| here and introduce the projected transverse-spin-overlap dis- 
tribution function -Prl'Zdiag) defined in terms of the diagonal overlap gaiag which is the trace of the tensor overlap 
Qfiv s, 

1 ^ (1) (2) (1) (2) 

tJ,=x.y i=l 

The distribution function PxC'Zdiag) is symmetric with respect to gdiag = 0. In the high-temperature phase, each g^jy 
(/i, v — x,y) is expected to be Gaussian-distributed around g^jy = in the L —t oo limit, and so is Qdiag- 

f21 

Let us hypothesize here that there exists a transverse spm-glass ordered state characterized by a nonzero , or 
by a nonzero EA transverse SG order parameter > 0. Reflecting the fact that gdiag transforms nontrivially under 
independent 0(2) rotations around the z-axis on the two replicas, which are the symmetries relevant to the transverse 
spin components in the presence of magnetic fields, even a self-overlap contributes nontrivial weights to Pxl'Zdiag) 
other than at ±gf^. In the L ^ oo limit, the self-overlap part of Pxl^diag) should be given by 

PT(gdiag) = Uiq^i,^) + ^ / =. (22) 

The derivation of Eq.(p2[) has been given in Ref. in the context of the XY SG. If the transverse SG ordered state 
accompanies RSB, the associated nontrivial contribution would be added to the one given by Eg. (p2|) . In any case, an 
important observation here is that, as long as the ordered state possesses a finite transverse SG LRO, the diverging 
peak should arise in PT('7diag) at (?diag = ±9t^ as illustrated in Fig.^, irrespective to the occurrence of the RSB. 







'^T(^diag) 












.qEA 

T 


gEA 
T 



FIG. 2. Sketch of the form of the transverse diagonal-spin-overlap distribution function PT(gdiag) in the thermodynamic 
limit, expected when there exists a finite transverse SG long-range order with a nonzero > 0. 
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IV. MONTE CARLO RESULTS 



This section is the core part of the present paper. Here, we present our MC results on the 3D ± J Heisenberg SG 
in magnetic fields. 



A. Chirality-related quantities 

First, we begin with the chirality-related quantities. In Fig.^ we show the temperature and size dependence of the 
mean local amplitude of the chirality for various fields. As can clearly be seen from FigH(a), extrapolation of xi^) to 
T = gives non-zero values as long as the applied field intensity is not too large, i.e., x{T ^ 0) ~ 0.294, 0.295, 0.308, 
0.313, 0.260, and 0.100 for H/J = 0, 0.1, 0.5, 2.0, 3.0, and 5.0, respectively. This indicates that the spin ordering 
of the 3D Heisenberg SG is certainly noncoplanar, which guarantees that the system sustains the nontrivial chirality. 
Meanwhile, a direct inspection of the spin pattern suggests that such noncoplanar spin configurations realized at low 
temperature in zero and weak fields is rather close to the coplanar one. Indeed, for completely random configurations 
of Heisenberg spins, x should take a value ■\/2/3 ~ 0.4714- • • a value considerably larger than the extrapolated 
x{T —^ 0) values. This again suggests that the noncoplanar configuration realized in zero and weak fields is close 
to the coplanar one. Interestingly, our MC data indicate that, in the weak field regime, x slightly increases with 
increasing magnetic field at fixed temperatures. This observation could be understood if one notes that the zero-field 
noncoplanar spin configuration is close to the coplanar one, and that the application of a magnetic field to such nearly 
coplanar spin configuration tends to "rise up" the spins from this plane with keeping the plane orthogonal to the 
applied field. This gives rise to more "three-dimensional" local spin structures with larger x- Of course, when the 
field is further increased, x eventually decreases simply because strong enough fields force spins to align along the 
field. In Fig.^(b), we show the size dependence of x for the field H/J = 0.1. As can be seen from the figure, there is 
very little size dependence in x- 
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0.3 
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0.1 
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H/J=0.5 


--■ 




H/J =1.0 


■o 




H/J=3.0 


-• 
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-A - - ' 



0.4 



0.3 
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1^ 0.2 



0.1 




H/J=0.l 







0.1 



0.2 



0.3 



0.4 







0.1 



0.2 



0.3 



0.4 



T/J 



T/J 



FIG. 3. Temperature dependence of the mean local amplitude of the chirality for various magnetic fields. The lattice size is 
L — 16 for H/J = 0, 0.1 and 0.5, and is L = 6 for other field values. For the case of H/J = 0.1, the size dependence of x is 
shown in Fig.(b). 



In Fig.0, we show the chiral-glass order parameter q^'' for the fields (a) H/J 
fields, increases rather sharply at lower temperatures. 



0.1, and (b) H/J = 0.5. For both 
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0.004 



0.003 



0.002 



0.001 



(a) H/J=0. 1 



L=6 ' — - 

L=8 • — 

L=10 ' — m- 

L=12 ' — e- 

L=16 ' — •- 



-m — e- 



0.1 



0.2 0.3 
T/J 



0.4 



0.004 



0.003 - 



0.002 



0.001 



(b) H/J=0.5 



L=6 

L=8 i — B- 

L=10 ' — ^- 

L=12 ' — 

L=16 ' — ^ 



0.1 



0.2 0.3 

T/J 



0.4 



FIG. 4. Temperature and size dependence of the chiral-glass order parameter for the field (a) H/J = 0.1, and (b) H/J — 0.5. 



In Figs||(a) and (b), we show the Binder ratio of the chirahty for the fields (a) H/J = 0.1, and (b) H/J = 0.5. 
As can be seen from the figures, exhibits a negative dip which, with increasing L, tends to deepen and shift toward 
lower temperature. Furthermore, g^ of various L cross at a temperature slightly above the dip temperature Tdip on 
negative side of g^, eventually merging at temperatures lower than T^ip. The observed behavior of g~^ is similar to 
the one observed in zero field [^. As argued in Ref. ||], the persistence of a negative dip and the crossing occurring 
at < are strongly suggestive of the occurrence of a finite-temperature transition where g^(TcQ) takes a negative 
value in the L oo limit. 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

T/J T/J 

FIG. 5. Temperature and size dependence of the Binder ratio of the chirahty for the fields (a) H/J = 0.1, and (b) H/J — 0.5. 



In Fig.^, we plot the negative-dip temperature T^ip{L) versus 1/L for the fields H/J = 0.1 and H/J = 0.5. For 
both fields, the data lie on a straight line fairly well. The linear extrapolation to 1/L = 0, as shown by the solid 
lines in the figure, gives our first estimates of the bulk chiral-glass transition temperature, i.e., Tqg/J — 0.23 for 
H/J = 0.1 and Tcg/J - 0.25 for H/J = 0.5. More precisely, rdip(L) should scale with L'^^" where f is the 
chiral-glass correlation-length exponent. As shown below, our estimate of ~ 1.3 comes close to unity, more or less 
justifying the linear extrapolation employed here. Indeed, extrapolation with respect to L^^^^-^, shown by the dashed 
curve in Fig.||, yields Tcg/J - 0.21 for H/J = 0.1 and Tcg/J - 0.23 for H/J = 0.5. 
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0.1 



H/J=0.5 ' — 



I ' ' ' 1 

0.05 0.1 0.15 0.2 

l/L 

FIG. 6. The dip temperature of Qx plotted versus l/L for the fields H/J = 0.1 and 0.5. L = oo extrapolation is performed, 
either based on the l/L fit (solid line), or on the l/L^/^-^ fit (dashed curve). The extrapolated values, Tdip(oo), give estimates 
of the bulk chiral-glass transition temperature Tea- 

As shall be argued below, we attribute several unusual features of g^, e.g., the growing negative dip and the crossing 
occurring at g^ < 0, to the possible one-step-like peculiar RSB in the chiral-glass ordered state. In systems exhibiting 
the one-step RSB, e.g., the mean-field three-state Potts glass, the Binder ratio is known to behave as illustrated 
in Fig.^ with a negative dip and the crossing occurring on the negative side |2^. Indeed, such a behavior is not 
dissimilar to the one we have observed in Fig.|5|. 



g ' 




1 




>N finite system 












FIG. 7. Sketch of the behavior of the Binder ratio expected in a system exhibiting a one-step RSB transition. The solid 
curve represents the behavior of an infinite system, while the dashed curve represents that of a finite system. 

An independent estimate of Tqq can be obtained from the equilibrium dynamics of the model. Thus, we also 
calculate the chirality autocorrelation function C^{t) defined by Eq.(p^). To check the possible size dependence, we 
show in Fig.|| the time dependence of C^{t) for the field H/J = 0.5 on a log-log plot, computed for (a) L — 16, and 
for (b) L = 20. As shown in the figures, C^it) shows either a downward curvature characteristic of the disordered 
phase, or an upward curvature characteristic of the long-range ordered phase, depending on whether the temperature 
is higher or lower than T/J ~ 0.23. Just at T/J ~ 0.23, the linear behavior corresponding to the power-law decay is 
observed. Hence, our data indicates that the chiral-glass transition takes place at Tqq/ J = 0.23(2), in agreement with 
our above estimate based on g^. From the slope of the data at T = Tqq, the exponent A characterizing the power-law 
decay of C^{t) « t^^ is estimated to be A = 0.13(2). We note that both our data of L = 16 shown in Fig.^(a) and of 
L = 20 shown in Fig.||(b) give almost the same estimates of Tqq, and of A, even though the L = 16 and L = 20 data 
themselves do not completely overlap, particularly below T/ J ^ 0.18. Anyway, our observation that C^(t) exhibits 
an upward curvature below Tcg, tending to a nonzero value corresponding to the static chiral EA parameter Qqq, 
indicates that the chiral-glass ordered state is "rigid" with a nonzero long-range order. The same analysis applied to 
the H/J = 0.1 case yields Tcg/J = 0.21(2) and A = 0.17(2). 
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FIG. 8. Temporal decay of the equilibrium chiral autocorrelation function at various temperatures for the field H / J = 0.5 
for the sizes (a) L = 16, and (b) L = 20. Temperatures correspond to T/J = 0.17, 0.18, 0.19, 0.20, 0.21, 0.22, 0.23, 0.24, 0.25, 
0.26 from top to bottom. Straight lines of power-decay fit are shown in both figures (a) and (b) at T/J = 0.23. 



In Fig.^ we show the chiral-overlap distribution function Px{Qx 
0.16, well below TcG 

grow and sharpen with increasing L, there appears a "central peak' 



for the field H/J — 0.5 at a temperature T/J = 

which 



In addition to the standard "side-peaks" corresponding to the EA order parameter ^Qqq , 

Stt Qx — for larger L, which also grows and 
sharpens with increasing L. The shape of the calculated P^i^x) very much similar to the one obtained in Ref. fqlin 
zero field, but is quite different from those observed in the standard Ising-likc models such as the 3D EA model ||27|] 
or the mean- field SK model |]l2|. As argued in Ref. in case of zero field, such peculiar features of PxiQx) are likely 
to be related to the one-step-like RSB. The existence of a negative dip in the Binder ratio gx and the absence of the 
standard type of crossing of gx at gx > are also consistent with the occurrence of such a one-step- like RSB [H,p3) . 
We note that our data of Px{qx) are also compatible with the existence of a continuous plateau between [— 9cq, 9cg] 
in addition to the delta-function peaks. 
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FIG. 9. The size dependence of the chiral-overlap distribution funcfi 
well below the chiral-glass transition temperature. Tea/ J — 0.23. 
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ion for the field H/J = 0.5 at a temperature T/J = 0.16, 



In FigJl^, we show the the temperature and size dependence of the Gx and Ax parameters for the field H/J ~ 0.5. 
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Although error bars of the data arc rather large here, the crossing occurs at temperatures somewhat higher than 
Tcg/J — 0.23 in both figures, while the crossing temperatures of neighboring sizes {e.g. L = 6 and L = 8 etc) 
gradually shift towards Tcg/J — 0.23 for larger L. The data are consistent with our estimate of Tqg above based on 
the Binder ratio and the autocorrelation. 
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B. Spin- related quantities 



In this subsection, we present our MC results of the spin-related quantities. In Figs.|Tl| and |T^, we show the spin 
Binder ratios for the longitudinal and transverse components, respectively, for the fields (a) H/J = 0.1, and (b) 
H/ J — 0.5. For both fields, the longitudinal Binder ratio gi^ increases monotonically toward unity with increasing 
L at all temperatures studied: See Fig.|ll|. This observation refiects the fact that the longitudinal component of the 
spin exhibits a net magnetization induced by applied fields at any finite temperatures. By contrast, the Binder ratio 
of the transverse component of the spin g^ decreases toward zero with increasing L, without a negative dip nor a 
crossing: See Fig. 12, This suggests that the transverse component of spin remains disordered even below Tqq. 
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FIG. 11. Temperature and size dependence of the Binder ratio of the longitudinal-component of the spin for the fields (a) 
H/J — 0.1, and (b) H/J = 0.5. For the field H/J = 0.5, magnified figure is shown in the inset. 
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FIG. 12. Temperature and size dependence of the Binder ratio of the transverse component of the spin for the fields (a) 
H/J = 0.1, and (b) H/J = 0.5. 



In Fig.[l^, we show the connected Binder ratio of the longitudinal spin component for the fields (a) H/J = 0.1, 
and (b) H/J = 0.5. Again, any anomalous behavior is not appreciable, no crossing nor extremum. Instead, 
monotonously approaches zero with increasing L, staying negative at any temperature. (Strictly speaking, the data 
of L = 6 and L = 8 for H/J ~ 0.1, exhibits a crossing-like behavior around T/J 0.2, but this is limited to these 
smaller lattices.) 
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In Fig. 14, we show the diagonal transverse-spin-overlap distribution function PT('i'diag) for the field H/J = 0.1 
at a temperature T/J = 0.18, well below the chiral-glass transition temperature Tqq = 0.21(2). The calculated 
^T(9diag) exhibits a symmetric "shoulder" at some nonzero value of qdiag, but as shown in the inset, this "shoulder" 
gets suppressed with increasing L, not showing a divergent behavior. Such suppression of the shoulder indicates that 
the chiral-glass ordered state does not accompany the standard transverse SG order, at least up to temperatures 
« (2/3)TcG- For H/J = 0.5, we have also observed similar suppression of the shoulder up to temperatures as low 
as around « (2/3)Tcg- Hence, we conclude that the chiral-glass ordered state does not accompany the standard 
transverse SG order, at least just below the chiral-glass transition point. Strictly speaking, the observed suppression 
of the shoulder is still not inconsistent with the Kosterlitz-Thouless(KT)-like critical SG ordered state. However, we 
note that such a critical SG ordered state appearing at T < Tqq is not supported by our data of shown in FigJl2. 









^diag 

FIG. 14. The diagonal transverse-spin-overlap distribution function -Pr(qdiag) for the field H/J = 0.1 at a temperature 
T/J = 0.18, well below the chiral-glass transition temperature Tcg/ J — 0.21. A magnified view of the shoulder part, indicated 
by the dashed circle in the main panel, is shown in the inset. 
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C. Critical properties of the chiral-glass transition 



In this subsection, we determine static and dynamical critical exponents associated with the chiral-glass transition. 
The analysis here is made for the two particular field values, H/J = 0.1 and H/J = 0.5, where most extensive 
simulations have been performed. In the analysis below, we fix Tcg to be Tqg/J = 0.21 {H/J = 0.1) and Tea/ J = 
0.23 {H/J = 0.5), as determined above. 

We estimate first the chiral-glass susceptibility exponent 7cg from the asymptotic slope of the log-log plot of the 
reduced chiral-glass susceptibility Xx versus the reduced temperature t = {T — Tcg) /Tcg- An example is given in 
Fig.[l5|(a) for the case of H/J — 0.5, where an asymptotic slope 7cg = 2.0(2) is obtained. 

We then estimate the chiral-glass critical-point-decay exponent TycG from the L-dependence of the chiral-glass order 
parameter q^'' at TcG, according to the relation q'^^ « j^-ii+ricG) ^ An example for the H/J = 0.5 case is shown in 
Fig.[l5|(b), where we plot q^ ^ at T/J ~ 0.23 ~ Tqg/J versus i on a log-log plot. As can be seen from the figure, the 
data lie on a straight line fairly well. From its slope ~ 1.5, the exponent rjcG is estimated to be rjcG = 0.5(3). 

The rest of the static exponents, acGj PcG^ and vcG^ can be estimated from 7cg and 77CG by using the standard 
scaling and hyperscaling relations as acG — ~l-9(5), PcG — 0.9(4) and vcG — 1.3(3). 
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FIG. 15. (a) Log-log plot of reduced chiral-glass susceptibility Xx versus the reduced temperature for the field H/J — 0.5. 
Its slope ~ 2.0 determines the chiral-glass susceptibility exponent 7cg = 2.0(2). The transition temperature is assumed here to 
be Tcg /J = 0.23. (b) Log-log plot of g^^' versus L for the field H/J = 0.5 at T/J = 0.23 ^ Tea/ J- Its slope ^ 1.5 determines 
the chiral-glass critical-point-decay exponent to be ?7cg = 0.5(3). 



The dynamical exponent zcg can be estimated from the exponent A, via the relation A — /ScgZ-^cg'^cg- From our 
above estimate, A = 0.13(2), we get zcg = 5.3(6). 

The same procedure is repeated for the case of H/J — 0.1. We then get ucG = 1-3(2), 77CG ~ 0.6(3), zcG = 4.9(6). 
These estimates for H/J = 0.1 agree within errors with the corresponding estimates for H/J — 0.5. Our estimates of 
the chiral-glass exponents are summarized in Table and are compared with the corresponding zero-field exponents 
reported in Ref. The finite-field exponents turn out to agree within errors with the corresonding zero-field 

exponents, suggesting that the zero-field and finite-field chiral-glass transitions lie in a common universality class. 
We note that this observation is consistent with the chirality scenario of Refs. 10. In Table ||, we also show 
the SG exponents of the 3D Ising EA model |27 2^ together with typical experimental values (in zero field) of real 
Heisenberg-like SG magnet AgMn ||2^. The critical properties of the chiral-glass transition differ clearly from those 
of the 3D Ising EA SG. By contrast, the chiral-glass exponents are close to the experimental exponent values for 
canonical SG AgMn, giving further support to the spin-chirality decoupling-recoupling scenario. 
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TABLE II. List of various critical exponents for the chiral-elass transition in zero- and in finite fields, compared with the 
corresponding spin-glass exponents of the 3D Ising EA model p7|j2q ] and of real Heisenberg-like SG magnet AgMn determined 
experimentally psl]. 





a 


/3 


7 


V 


T] 


z 


chiral-glass: H/ J = 0.1 


2.0 


- 1.0 


- 1.8 


- 1.3 


- 0.6 


-4.9 


chiral-glass: H/ J = 0.5 


2.0 


- 0.9 


- 2.0 


- 1.3 


- 0.5 


-5.3 


chiral-glass: H/ J = 0^^ 


1.7 


- 1.1 


- 1.5 


- 1.2 


- 0.8 


- 4.7 


3D Ising EA27'28) 


- -1.9 


0.55-0.67 


4.0-4.7 


1.7-2.0 


-0.37—0.26 


6.0-7.0 


AgMn^^^ 


-2.2—1.9 


0.9-1.0 


2.1-2.2 


1.3-1.5 


0.3-0.4 


- 5.4 
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As a consistency check of our estimates of exponents and Tcg values, we have also done the following: We use 
the ucG value determined above, ucg ~ 1-3, and extrapolate the dip temperature of g^, T^ip{L), to L = oo (see 
the dashed lines of Fig.^) . As mentioned, such an extrapolation yields the bulk chiral-glass transition temperature, 
Tcg = 0.21(2) {H/J = 0.1) and Tcg = 0.23(2) {H/J = 0.5). These estimates of Tec agree with those obtained from 
the chiral autocorrelation and employed in our scaling analysis. This guarantees that our analysis of exponents and 
Tcg is self consistent. 

In Fig.^, we show the the standard finite-size scaling plot for the chiral-glass order parameter q'^^ based on the 
relation, 

«i-(i+')cG)j(2,t-i/'^cG) ^ (23) 

where the Tcg, ''7CG and i^cG values are set to the best values determined above. As can be seen from the figures, 
reasonable data collapsing are obtained, at least for larger lattices. At the same time, however, one sees that there 
exists a systematic deviation from the scaling for smaller lattices, particularly in the case of H/J ==0.1. Such a 
deviation observed for smaller lattices suggests the existence of a significant finite-size correction. 
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FIG. 16. Finite-size plot of the chiral-glass order parameter for the fields (a) H/J — 0.1, and (b) H/J = 0.5 
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The existence of such significant finite-size effects has also been suggested from the behavior of the chiral-overlap 
distribution function Pxilx) ^'^'^ '^^ chiral Binder ratio g^^. In a truly asymptotic critical regime, Px{qx) itself 
should scale at T = Tcg with tuning one exponent 77cG- However, in the range of sizes studied here L < 16, we 
cannot observe such a full scaling of P-^{q^). Such a lack of complete scaling of Pxilx) gives rise to certain degrees 
of uncertainty in our estimate of rycc '■ Namely, if one tries to scale the width of the distribution such as its second 
moment qx \ it yields ?7cG ~ 0.5 as given above (see Fig. |l^(b)), while if one tries to scale the height of P^iqx), it 
instead yields 77CG ^ 0.4, which is somewhat smaller than the above estimate, though still lying within the quoted 
error bar. Lack of a complete scaling in P^iqx) is also reflected in the behavior of g^, which does not show a unique 
crossing at T = Tec within the range of sizes studied: Instead, as shown in Fig.^, the crossing occurs on the negative 
side of g^ considerably above T = Tcg, while the crossing points tend to come down toward Tcg as L increases. 

Concerning the transverse spin order, from the behaviors of the Binder ratio and of the diagonal transverse-spin- 
overlap distribution function, we have already found a strong numerical evidence that the chiral-glass transition does 
not accompany the transverse SG order, at least just below Tcg- In other words, the transverse component of the spin 
orders only at zero temperature, or else, if it orders at a finite temperature, the associated transverse SG transition 
temperature Tgc is significantly lower than Tec, say, below « (2/3)Tcg- We warn the reader here that, so long as 
one looks at the SG correlation or the SG order parameter, a rather careful analysis is required to really see such a 
behavior. As an example, we show in Fig.p^the standard finite-size scaling plots of the transverse SG order parameter 

(2) 

Qj for H/J ~ 0.1; (a) the one assuming Tsc = 0, and (b) the other assuming Tsc = 0.21J(= Tcg)- Similar plots are 
given in Fig.|l| for the field H/J = 0.5 with assuming Tsc = 0, assuming (a) Tsc = and (b) Tsc = 0.23J(= Tcc)- 
At a look, both fits seem equally acceptable without appreciable difference if the exponents are adjusted in appropriate 
ways. Then, one may wonder if the transverse spin might order simultaneously with the chirality, with the associated 
SG exponents i^sG — 1.1(~ J^cc) and r]sG — —0-25: See Fig.^ We believe, however, this not to be the case due to 
the following reasons. 
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FIG. 17. Finite-size scaling plot of the transverse spin-glass order parameter for the field H/J — 0.1, assuming (a) Tsa = 0, 
and (b)TsG = 0.21J(=rcG). 
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FIG. 18. Finite-size scaling plot of the transverse spin-glass order parameter for the field H/J — 0.5, assuming (a) Tsg = 0, 
and (b) Tsg = 0.23 J(= Tog). 



First, as shown above, such simultaneous chirahty and transverse-spin ordering contradicts with our result of the 
Binder ratio and the diagonal transversc-spin-overlap distribution function -Prl'Jdiag)- Second, a closer inspection 
of the data reveals that, at and below T — Tqq, there exists an important difference between the behaviors of the 
transverse spin q^^ and of the chirality q'^^ . 

In Fig.|l^, we show on a log- log plot the size dependence of the both order parameters, q!^'' and q^K at several 
temperatures at and below Tqg ■ As can be seen from Fig.|l^(a), below Tqg the chiral-glass order parameter exhibits 

(2) 

a clear upbending for larger L, indicating that qj, tends to a nonzero value in the thermodynamic limit. In sharp 

(2) (2) 

contrast to this, such an upbending is never seen in the transverse SG order parameter q^, . Instead, qj, shows a 
slight downbedning behavior at T = Tcg, which gradually shifts to the near linear behavior at lower temperatures. 

(2) 

The observed behavior of is consistent with either, (a) the onset of the Kosterlitz-Thouless(KT)-like transition at 
a finite temperature below which the spin-glass correlations decay algebraically with a power-law or, (b) the gradual 
growth of the transverse SG correlation length ^ which exceeds the investigated system size L — 16 around a certain 
nonzero temperature close to Tcg- In the former case, there should exist a well-defined finite SG transition temperature 
with the critical SG ordered state, while, in the latter case, there need not be a thermodynamic SG transition at a 
finite temperature. Generally speaking, it is difficult to discriminate between the above two possibilities only from 

(2) 

the data of finite sizes with L < ^. 

Nevertheless, we believe we can at least exclude here the possibility that the KT-like transverse SG transition 
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occurs simultaneously with the chiral-glass transition at T = TqGi accompanied by the critical SG ordered state at 



T < TcG- First, we note that such a critical SG ordered state is not supported by our data of gT of Figjl2[ Second, 
the transverse spm-glass correlation-length exponent estimated in FigsjT7|(b) and ^(b) assuming the simultaneous 
spin and chiral transition, v ~ 1.1, is far from from the lower-critical-dimension (LCD) value, v = oo, generically 
expected for such a KT-like transition. In so far as one insists that the transverse SG order occurs simultaneously 
with the chiral-glass order, our numerical estimate of the transverse SG correlation-length exponent is not compatible 

(2) 

with the LCD value v — oo, which is now hard to reconcile with the KT-like behavior observed in grp at T < TcG- 
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FIG. 19. Size dependence ol the order parameter for the field H/J — 0.1, at several temperatures at and below Tcg'- 



(a) the transverse-spin-glass order parameter q!^' , and (b) the chiral-glass order parameter q'^\ The chiral-glass transition 
temperature at this field is Tcg/ J — 0.21. To emphasize the deviation from the linearity, lines connecting the two small-size 
data L — 6 and 8 are drawn at each temperature. 



In fact, as recently argued in Ref. |g6[ for the case of the XYSG, the chirality scenario gives the possible cause 
why simultaneous spin and chiral orderings are apparently observed in the SG order parameter or the SG correlation 
function. This would closely be related to the length and time scales of the measurements. Here, one should be 
aware of the fact that the spin-chirality decoupling is a long-scale phenomenon: At short scale, the chirality is never 
independent of the spin by its definition, roughly being its squared (x ~ S^Sj^ S^) as expected from the naive power 
counting. Hence, the behavior of the spin-correlation related quantities, including the SG order parameter which is a 
summed correlation, might well reflect the critical singularity associated with the chirality i.e., the one of the chiral- 
glass transition, up to certain length and time scale. In such a scenario, apparent (not true) transverse "spin-glass 
exponents" expected would be v^q ~ i^cG ^ 1-3 and rj'gQ ^ —0.25, the latter being derived from the short-scale 
relation, 1 -I- 77CG ~ 2(1 -I- Tj'gn). Note that these values are not very far from the ones we get from the finite-size scaling 
analysis of Figs.|l^(b) and |8|(b), assuming the simultaneous occurrence of the spin and chiral transition. However, 
we stress again that such a disguised criticality in the spin sector is only a short-scale phenomenon, not a true critical 
one. 



D. Phase diagram 



By collecting our estimates of the Tcg values for various field values, as obtained by the extrapolation of TdiJJj) to 
L — 00, we construct a phase diagram in the temperature vs. magnetic field plane. The result is shown in FigJ2n. We 
have used here the zero-field estimate of Ref. |30|, Tcg/ J ~ 0.21(2). Error bars are estimated here from the differences 
between the extrapolated Tcg values via the l/L and l/L^^^-^ fits. As is evident from Fig.^ the chiral-glass state 
remains quite robust against magnetic fields. This is most evident in Fig.|2^(b) where we draw the same phase diagram 
on a plot where both the temperature and the magnetic-field axes have common energy scale. Indeed, Tgg{H) is not 
much reduced from the zero-field value even at a field as large as ten times of Tcg(O). At lower fields, the chiral-glass 
transition line is almost orthogonal to the H = axis, consistent with the behavior Eq.(^ derived from the chirality 
scenario. Our data are even not inconsistent with the coefficient c in Eq.(||) being slightly negative so that Tcg{H) 
initially increases slightly with H, though it is difficult to draw a definite conclusion due to the scatter of our estimate 
of Tcg{H). If one remembers here our MC observation that the application of a weak magnetic field tends to increase 
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the mean local amplitude of the chirahty, x, from its zero-field value, such an initial increase of Tcq{H) seems not 
totally unlikely. 
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FIG. 20. H-T phase diagram of the 3D ± J Heisenberg SG determined by the present simulation. Note that the energy scales 
of the H and of the T axes are mutually different in Fig. (a), while they are taken to be common in Fig.(b). 



V. SUMMARY AND DISCUSSION 



In summary, we have performed large-scale equilibrium Monte Carlo simulations on the 3D isotropic Heisenberg SG 
in finite magnetic fields. We have confirmed that our MC results are consistent with the chirality scenario of Ref. [|. 
Among other things, we have verified the occurrence of a finite-temperature chiral-glass transition in applied fields, 
essentially of the same character as the zero-field one. The chiral-glass ordered state exhibits a onc-stcp-like peculiar 
RSB, while it does not accompany the transverse SG order, at least up to temperatures around « (2/3)Tcg- The 
criticality of finite-field chiral-glass transitions seems to be common with that of the zero-field one, which, however, 
clearly differs from the criticality of the standard 3D Ising EA model. Meanwhile, the chiral-glass exponents turn out 
to be close to the experimental exponents determined for canonical SG such as AgMn. We have also constructed a 
magnetic phase diagram of the 3D Heisenberg SG model. The chiral-glass transition line in the H-T plane is found 
to be almost vertical to the temperature axis, up to rather high fields of order H ~ IOTcg(O), indicating that the 
chiral-glass ordered state is quite robust against magnetic fields. This somewhat surprising property probably arises 
from the fact that the magnetic field couples in the Hamiltonian directly fo the spin, not to the chirality, and the 
effective coupling between the field and the chirality is rather weak. The chiral-glass transition line has a character 
of the Gabay- Toulouse line of the mean-field model, yet its physical origin being entirely different. 
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FIG. 21. Experimental phase diagram of AuFe as determined by Campbell et al by means of torque measurements [Ref. 
[0]. For comparison, we also show the present numerical result of the magnetic phase diagram of the 3D ±J Heisenberg SG 
model. The way how we scale the units of magnetic field and temperature in plotting the experimental data is explained in 
the text. 



It is not immediately possible to make a direct comparison of our results with experiments. This is mainly because 
the random magnetic anisotropy, which inevitably exists in real SG materials, is not introduced in our present model. 
Furthermore, in real SG magnets, spins do not necessarily sit on a simple-cubic lattice, nor interact with other spins 
via the nearest-neighbor ±J coupling, etc. In spite of these obvious limitations, it might be interesting to try to 
compare our present magnetic phase diagram with the experimental one for Heisenberg-like SG magnets. Chirality 
scenario claims that, in the high-field region where the anisotropy is negligible relative to the applied magnetic field, 
the SG transition line should essentially be given by the chiral-glass transition line of the fully isotropic system. If so, 
our present result entails that the SG transition line of real Heisenberg-like SG should be almost vertical against the 
temperature axis in the high-field regime where the magnetic field overwhelms the random magnetic anisotropy. In 
Fig.^(a), we reproduce the experimental H-T phase diagram of canonical SG AuFe from Ref. [|ll|. In the same figure, 
we also show our present result of the chiral-glass transition line, scaled in the following way. We try to mimic the real 
system by the classical Heisenberg Hamiltonian with an effective coupling J and an effective magnetic field , which 
is defined in terms of Eq.(^). First, we estimate the zero-field transition temperature of the hypothetical isotropic 
system to be Tg k, lOK, by extrapolating the high-field GT-like transition line of AuFe to H ~ Q. Then, with the 
knowledge of our present estimate of Tcg ~ 0.2 J, we estimate the relevant J roughly to be 50K. The field intensity H 
is then translated into the field intensity in the standard unit H* by the relation H = pcffH* , peff being the efi^ective 
Bohr number: In case of AuFe, Peff was experimentally estimated to be 4.55/xb, where fj,B is the Bohr magneton 
Thus, our Fig.^ suggests that the SG phase boundary of AuFe might stay nearly vertical up to the field as high 
as H lOrsG(O) ^ 40[T]. Of course, considering the difference in microscopic details between the present model 
and real AuFe, one cannot expect a truly quantitative correspondence here. Anyway, further high-field experiments 
on AuFe and other Heisenberg-like SG magnets might be worthwhile to determine the SG phase boundary in the 
high-field regime. 

In order to make further comparison with the experimental phase diagram in the low-field regime, it is essential to 
examine the efi^ects of random magnetic anisotropy inherent to real SG materials. Indeed, in the low-field regime where 
the applied field intensity is comparable to or weaker than the random magnetic anisotropy, the chirality scenario 
predicts the appearance of a singular crossover line which has some character of the AT-line of the mean-field model 
1^,^. In order to make further insight into the spin-glass and the chiral-glass orderings in magnetic fields and to 
check further the validity of the chirality scenario, it would be interesting to make similar finite-field simulations for 
the anisotropic 3D Heisenberg SG model. 
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